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ABSTRACT 

We present DLSCL J0916. 2+2951 (z=0.53), a newly discovered major cluster merger in which the 
coUisional cluster gas has become dissociated from the collisionless galaxies and dark matter. We 
identified the cluster using optical and weak lensing observations as part of the Deep Lens Survey. 
Our follow-up observations with Keck, Subaru, Hubble Space Telescope, and Chandra show that the 
cluster is a dissociative merger and constrain the dark matter self-interaction cross-section ctdm^dm ^ 
7cm^g~^. The system is observed at least 0.7 ± 0.2 Gyr since first pass-through, thus providing a 
picture of cluster mergers 2-5 times further progressed than similar systems observed to date. This 
improved temporal leverage has implications for our understanding of merging clusters and their 
impact on galaxy evolution. 

Subject headings: cosmic background radiation — dark matter — galaxies: clusters: individual (DLSCL 
J0916. 2-1-2951) — gravitational lensing: weak — X-rays: galaxies: clusters 



1. INTRODUCTION 

During a major galaxy cluster merger the collision- 
less dark matter (DM) and galaxies become disso- 
ciated from the cluster gas (the principal baryonic 
component), a process that is observable for ^2 Gyr 
(|Ricker fc SarazinI [200l . These mergers have become 
important astrophysical probes providing a diagnostic 
of large scale structure (LSS) formation, e mpirical evi- 
dence favoring DM over modified gravity (jClowe et al.l 
[20061) . constraints on the DM self-interaction cross- 
section ([Markevitch et aLll2004D . and c onstraints on th e 
large scale matter-antimatter ratio (jSteigmanl I2008D . 
They are als o a suspected source of extremely energetic 
cosmic rays (jvan Weeren et al1l2010( ) . A special subclass 
of mergers are probes of all the aforementioned science. 
Such mergers: 1) are between two subclusters of com- 
parable mass, 2) have a small impact parameter, 3) are 
observed during the short period when the cluster gas 
is significantly offset from the galaxies and DVI, and 
4) occur mostly transverse to the line-of-sight so that 
the apparent angular separation of the cluster gas from 
the galaxies and DM is maximized. We term these sys- 
tems dissociative mergers. Given these requirements it is 
not surprising that o nly six such system s have been con 
firmed : IE 0657-56 (IClowe et al.l [2006l). MACS J0025.4 



1222 (iBrada c et al.l 120081). Abeh 520 (iMahdavi et al 



2007t [Okabe fc Umetsd 120081). Abell 2 744 (iMertenejal 



20Tll). AbeU 2 163 fOkabe et al.l[2MTI). and Abeh 1758 
(jRagozzine et a l. 2011). With such a small sample which 
only probes a narrow temporal window (^0.17-0.25 Gyr 
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post-collisior0), it is not surprising that we still do not 
completely understan d merging systems. For example, 
IMahdavi eFall (|2007f ) and Jee et al. (2012, ApJ, in press) 
find a "DM core" coincident with X-rays but not with any 
bright cluster galaxi es in the merging system A520, and 
iMerten et aH (|20lH ) find a "ghost" clump of gas lead- 
ing a DM clump in A2744; there is no accepted expla- 
nation for either of these observations. Simulations are 
needed to provide insight into these mysteries, but hydro- 
dynamic simulations are complex and many assumptions 
are re quired (see e.g. lSpringel fc Farraril2007l : lPoole et al.l 
|2006( ): therefore observations of new and different sys- 
tems will play a crucial role in building confidence in 
these simulations. 

We have identified a new dissociative merger, DLSCL 
J0916. 2-1-2951, that probes an unexplored area of merger 
phase-space. We originall y detected the cluste r in the 
Deep Lens Survey (DLS; IWittman "eFall I2OOI via its 
weak lensing (WL) shear signal. It consists of two main 
subclusters spectroscopically confirmed to be at the same 
redshift (0.53). This cluster was also observed in the 
Sunyaev-Zel'dovich Array Survey ([Muchovei et al.ll2011l) 
which provided evidence that the cluster gas is dissoci- 
ated from the bulk of the mass and galaxies. Follow-up 
optical observations with Subaru and HST enable higher 
resolution mass maps and follow-up X-ray observations 
with Chandra ACIS-I confirm that the majority of the 
gas is offset between the North and South subclusters, 
the signature of a dissociative merger (Figure [T]). 

In this letter we introduce DLSCL J0916.2-f 2951 and 
summarize our survey of its three dominant components 
(galaxies, DM, and gas) and the cluster's astrophysical 
implications. A more thorough exposition of the survey 
and analysis will be presented in Dawson et al. (in prepa- 
ration). Throughout this paper we assume f^A = 0.7, 
flm = 0.3, and H = 70kms-iMpc~^ 

^ We define the time of collision to be the time of the first peri- 
centric passage. 
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Fig. 1.— Merging cluster DLSCL J0916. 2+2951 and its three 
matter components. Overlaid on the HST color image of the galax- 
ies is the total mass distribution (blue) based on WL analysis of 
the HST images and the cluster gas distribution (red) based on 
Chandra X-ray observations. The bulk of the collisional gas is 
located between the two collisionless galaxy and mass concentra- 
tions, indicative of a dissociative merger. The image is 5' X 5' 
(~ 1.9 X 1.9Mpc2 at z = 0.53). 

2. OPTICAL ANALYSIS 

We obtained spectroscopic redshifts for 20 cluster 
members with Keck LRIS (2007 January 16) and 634 
unique spectroscopic redshifts (0 < z < 1.2) in a ~ 15' x 
15' area centered on the cluster with Keck DEIMOS 
(2011 March 2-3), including 132 members at the clus- 
ter redshift. We reduced the LRIS spectra using a 
scripted sequence of standard IRAF reduction tasks, 
and the DEIMOS spectra using a modified ver sion of 
the pEEP2 spec2d packag e (jDavis et al. 2003; Gal et al.l 
120081 : Eimux et al.ll2009t ). 

We use our full sample of 654 spectroscopic redshifts as 
well as photometric redshifts to identify potential line-of- 
sight structures which may confuse our results. We find 
no evidence for significant line-of-sight structure (Figure 

ED. 

We estimate each subcluster's redshift and velocity dis- 
persion (Table [ij using the biwe ight-statistic and bias- 
corrected 68% confidence limit () Beers et al.l I199O0 ap- 
plied to 100,000 bootstrap samples of each subcluster's 
spectroscopic redshifts. Our redshift estimates indicate 
a line-of-sight velocity difference of vios = 670^33" km 
s~^ between the North and South subclusters, using the 
galaxies within a 0.5 Mpc radius centered on the HSTWL 
mass peaks and within a velocity range of ±3000 km s~^ 
of 2=0.53 ('^ 3x the expected velocity dispersion); corre- 
sponding to 38 and 35 galaxies for the North and South 
subcusters, respectively. These results are robust against 
varying the velocity range ±1000 km s"""^ and using the 
Subaru WL or galaxy number density peaks as the aper- 
tures' centers, provided the aperture radius is <0.5Mpc: 
larger radii lead to significant subcluster membership 
confusion. Additionally, we report the velocity disper- 
sion mass es t imate s based on the scaling relation of 
lEvrard et all pOOSi ) in Tabled 

We obtained B, V, R, and z' photometric data (12, 
12, 18, and 12ksec, respectively) with Mosaic 1 on the 



KPNO 4-m May all telescope as part of the DLS. To 
improve the accuracy of our photometric redshifts we 
also observed the cluster in three medium- width optical 
bands {g,h, and i from the BATC filter set), bracket- 
ing the redshifted 4000 A feature, using the upgraded 
Mosaic 1.1 imager on the KPNO May all with exposure 
times of 6ks per filter (201 1 April 22-24). We esti- 
mate colors using ColorPro (Co e et al.l l2006l ) and red- 
shifts using BPZ (Bcnftez 2000). We replace the stan- 
dard templates with a s et "tweaked" in a method simi- 
lar to that described in lllbert et al.l (120061) . u sing spec- 
troscopic samples fro m SHELS dCeller et al] r2005') and 
the PRIMUS survey (|Coil et al.ll2010D which overlap the 
DLS. Figure [2] shows the density isopleths of galaxies 
with 0.43 < Zphot < 0.63 (roughly the cluster redshift 
i^'zphot)- This map agrees well with the distribution of 
spectroscopically confirmed cluster members. 



3. WEAK LENSING ANALYSIS 

To map the total mass distribution we use a version of 
the Fische r fc Tvso n (1997) method modified to include 
a novel tomographic signal-matched filter. The cluster's 
WL shear signal, 7, depends not only on the projected 
surface mass over-density of the cluster, AS, but on the 
relative distances of the observer, the mass, and the back- 
ground galaxies: 



7 



AS 



Ds{Zs) 



Zl 



where Scr is the critical surface density, T-l is the Heav- 
iside step function, and Di, Dg, & Di^ are the angular 
diameter distances to the lens, source, and between the 
lens and source, respectively. Since we do not have ex- 
act redshift measurements of the source galaxies we use 
each galaxy's photometric redshift probability distribu- 
tion function, p(z), to estimate a respective Scr for a 
given lens redshift, 

^crizi) ~ {^crizi)) = J T,cr{zi , Zs)p{Zs)dZs . 

In addition to weights based on shape measurement er- 
rors, we also weight the shear of each galaxy based on its 
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/ [T.cr{zi,Zs) - (Scr(zi))] p{Zs)dZs 



This method incr eases the signal-to-noise o f the mea- 
surement (see e.g. iHennawi" fc SDergelll2005l) . and more 
accurately accounts for the errors inherent in the pho- 
tometric redshift estimates, compared to single-point es- 
timates. We estimate uncertainties using 100 bootstrap 
resamplings. 

Encouraged by the DLS mass and galaxy maps we ob- 
tained higher-resolution ground and space based obser- 
vations. Subaru Suprime-Cam j'-band coverage of the 
cluster was provided by engineering-time observations of 
DLS Field 2 (2008 January 8). We use the Suprime- 
Cam data reductio n software SDFRED (jYagi et al.ll2002l : 
lOuchi et all }200l followed by SCAMP & SWARP 
( Bertin et all 120021: lBertinll2006l ) to refine the astrome- 
try and make the final mosaic. DLSCL J0916. 2+2951 
was also observed with HST ACS/WFC using F606W 
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Fig. 2. — Left: DLS composite BVR color image of DLSCL J0916. 2+2951 sliowing the galaxies of the two subclusters. The white contours 
represent the number density of galaxies with Zpij^t = 0.53 it 0.1, the cluster redshift io"zpi,ot ■ '^^^ contours begin at 200 galaxies Mpc~^ 
with increments of 50 galaxies Mpc~'^. The image field-of-view is the same as FigurefT] Right: Histogram of the 200 observed spectroscopic 
redshifts within the field of view of the left figure. The red portion is the subsample that passes the Zpi^^f^ — 0.53 it 0.1 criteria. The galaxies 
at 2 ^ 0.6 had equal probability of selection as the cluster members and show no sign of clustering. 

aru for the surrounding area. Using a tomography-based 
MCMC analysis we simultaneously fit NFW halos cen- 
ter ed on the North and Sou th HST WL peaks, and use 
the iGelman fc RubinI (|1992f) convergence test applied to 
eight independent chains. In order to reduce the number 
of free parameters we use the lDuffv et all ([200l empiri- 
cal relation between M200 and concentration. We present 
the most likely masses for each halo along with the bias- 
corrected 68% confidence limits in Table [T] We also 
compare the integrated projected surface mass density 
of the NFW hal os with the measured WL aperture mass 
(|Fahlman et al . 1994) of each subcluster and find general 
agreement within a radius of 0.5 Mpc of each subcluster. 

4. X-RAY ANALYSIS 

We acquired X-ray spectral-imaging of the cluster with 
40ks of Chandra ACIS-I time (GO-12800854), and re- 
duce it using CI AO version 4.2 and CALDB version 
4.4.1. We manually identify X-ray point sources and 
mask them before adaptively smoothing the diffuse emis- 
sion, we present the resulting map in Figure |31 We esti- 
mate source counts and their error using the dmextract 
function of CIAO. We use an 8' radius background region 
which encloses the subcluster regions and rests ~ 90% 
on ACIS-I3 (on which the subclusters are observed); yet 
excludes subcluster, chip gaps, and point sources. For 
the South and Central X-ray concentrations (120±17 and 
170±19 detected 0.5-2 keV photons, re spectively) we use 
the Xspec X-ray spectral fitting tool ([Arnaudl 119960 to 
fit a Mewe-Kaastra-Liedahl plus photoelectric absorp- 
tion model (fixed t o the Leiden/ Argentine/Bonn value; 
iKalberla et al.l26 0li) to the X-ray spectrum of each X-ray 
concentration. For the North concentration there are not 
enough detected X-ray photons (38±12) to fit a mean- 
ingful spectrum. We report the results of this analysis in 
Table [J 



and F814W fihers (GO-12377, PLW. Dawson) in a 2 x 1 
pointing mosaic that covers the subclusters (Figure [1]). 
The exposure times for F606W and F814W are 2520s 
and 4947s per pointing, respectively. We reduce this 
data followin g a rn ethod similar to that presented in 
iJee fc Tvsoi^ (|2009t ). We measure the PSF of both 



datase ts using the PCA method presented in iJee et al.l 

(IlooT). 

We perform our WL analysis independently on both 
the Subaru and HST F814W data. The Subaru data 
has 0.72" seeing and 49 WL-quality source galaxies 
(i.e. background galaxies with measured ellipticity error 
< 0.3) per arcmin^. For the mass map we use an apodiz- 
ing kernel radius of 0.5', which can be interpreted as the 
effective resolution of the WL mass map. We are able to 
cross-match most of the detected objects with the DLS 
and use the p{z)'s discussed in the previous section. 

Cross-matching is more problematic with the higher- 
resolution HST data, so we use a color-magnitude cut 
(F606W-F814W<0.8 and 24 < F814W < 28.5) to select 
source galaxies and exclude cluster red-sequence and 
bright foreground galaxies. For the WL analysis of the 
HST F814W image, which has a 0.1" PSF and 136 WL- 
quality source galaxies per arcmin^, we use an apodiz- 
ing kernel radius of 3.6". We find no significant spatial 
correlation between source density and subcluster posi- 
tion, suggesting that our source galaxy population is not 
significantly contaminated with cluster galaxies. We esti- 
mate the p(z) of our HST source galaxy sample by assum- 
ing the photometric redshift distribution of the lCoe et al.l 
((2OO6. ) HUDF catalog after applying our color-magnitude 
cut. Figure [3] shows general agreement between the Sub- 
aru and HST WL mass, and galaxy density maps. 

We construct a joint catalog from the HST and Sub- 
aru data, using the HST data where available and Sub- 
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Fig. 3. — Comparison of the Subaru i'-band ground-based (left) and HST space-based (right) WL mass signal-to-noise maps (color) of 
DLSCL J0916. 2+2951 with the X-ray distribution (bold black contours) and galaxy number density (white contours, same as Figure [2]l. 
The peak centers and corresponding one sigma errors are denoted by the gray cross-hairs. In both analyses there is general agreement 
between the location and relative magnitude of galaxies and WL yet the majority of the cluster gas is centered ~ 1.4' between the North 
and South subclusters in a local mass underdensity, providing evidence that the North and South subclusters have undergone the first 
pass-through of a major merger. The scale of each map is equivalent and the image field-of-view is the same as Figures [T] & [2] The map 
created from the joint Subaru/HST catalog looks nearly identical to the HST map, with only slight variations in the scale (see Table \Tl . 



5. CLUSTER MERGER SCENARIO 

The peak of the gas distribution (09^6^15! ± 
5.5!, 29?50.'59'.'±5.0'.') derived from X-rays is offset 1.4' ± 
0.49 from the North HST WL mass peak (09^6^10! ± 
30!, 29?52fl0'.'±30'.'), and 1.4'±0.14 from the South HST 
WL mass peak (09l'16^15! ± 8.0!, 29?49f34'.' ± 6.9'.'), and 
is located near a local minimum in the mass, suggesting 
that the subclusters have a small impact parameter and 
have experienced at least their first pass-through along a 
north- northwest merger direction (see Figure [J]). Further 
evidence for the merger scenario is provided by the mor- 
phology of the gas. S imulations (jSchindler fc Muelleii 
ll993HPoole et alJ 120061 ) predict that the gas morphol- 
ogy elongates transverse to the merger direction after 
pass-through for mergers with small impact parameters. 
The Central gas concentration appears to be oblate and 
roughly perpendicular to the axis connecting the North 
and South mass peaks. This is consistent with the in- 
terpretation that these two subclusters have experienced 
their first pass-through and that the merger axis being 
roughly in the plane of the sky. 

We gain further insight into the dynamics of the sys- 
tem by studying the spectroscopic redshifts of the sub- 
cluster galaxies. The observed line-of-sight velocity dif- 
ference, wios(iobsorvcd) — 670^33g km s""'^, between the 
North and South subclusters does not explicitly deter- 
mine the dynamics or orientation of the system since 
there is a degeneracy between the three-dimensional rel- 
ative peculiar velocity, I'reiativej and angle of the pair 
with respect to the plane of the sky, a. In order to con- 
strain this degeneracy we perform a Monte Carlo analy- 
sis of the system assuming the clusters to be uniform- 
density spheres with i?200 radii. We draw randomly 
from each subcluster's mass probability distribution as 



determined by the NFW MCMC analysis, the redshift 
distribution of the spectroscopic bootstrap analysis, and 
the centroid distribution as determined by the bootstrap 
sample of the HST WL mass map. From these parame- 
ters we calculate a conservative upper bound on the rel- 
ative velocity at the time-of-collision, WrGiative(icoiiision), 
by assuming that the two subclusters freely fell from 
an infinite separation and zero relative velocity yet be- 
gan with their current mass. To determine a conserva- 
tive lower bound we calculate Wroiativo(^coiiision) assum- 
ing the subclusters had reached their observed separa- 
tion with Wrciativc(iobscrvcd) = 0. Implicit in our anal- 
ysis is that the merger occurs with an impact parame- 
ter of zero. Including a moderate impact parameter of 
^0.1i?20Q shou ld have less than a 1% eff ect on the merger 
velocity (Mas tropietro fc Burkei^l2008| ) . As can be seen 
from Figure|4]the derived time-since-collision is relatively 
insensitive to moderate changes in ^relative (^collision 

). For 

a consistency check we compare our analytic estimate of 
the Bullet Cluster's t ime-since-collision with t he numer- 
ical model results of Spring el &: Farraii ([2007) and find 
agreement within 20% (compare the square and diamond 
points of Figure E]). While we attempt to be conserva- 
tive in our analysis, we disregard many complexities of 
the merger. Further insight will require deeper X-ray ob- 
servations and/or hydrodynamic numerical simulations. 

Using our analytic model we calculate that 
Wrciativc(icoiiision) = 17001280 s"^ gi'^eu the sound 



speed of the gas 
Mach number of 



600 km s ^, this corresponds to a 
3. With a projected separation 



of I.OI0.14MPC and estimated a = 341^° degrees we 



find a physical separation of 1.3_Q'^g Mpc between 
the North and South subclusters. This is the largest 
physical offset for a dissociative merger observed 
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TABLE 1 

Observed subcluster and X-ray concentration properties 



Subclustcr 


Redshift 


(km s-l) 


0"u M2OO 

(IO^Mq) 


WL M200 
(1O"M0) 


^^0.5-2kcV 

(10*^ergs-i) 


Tx 
(kcV) 


X-ray 
S/N 


Joint WL 
S/N 


North 
South 
Central 


53414+0 00065 


770111° 


3.7 ±2.3 
4.1 ± 1.6 


1 7-1-2.0 

^■'-0.72 
.+1.2 
■^•^-0.79 


0.63 
2.1 
2.8 


9 7-1-1.2 
^- '-0.7 
2 2+1'* 


3.2 
7.0 
9.1 


3.0 
6.7 

-3.3=^ 



The negative WL S/N indicates a projected surface mass local under-density. 



to date. This translates into a time-since-collision 
(^observed " ^collision) of O.TIq i Gyr. Give n the collision 
veloc i ties of the other dissociat i ve mergers (iBradac et alJ 
2001: iMarkevitch et all 12001 iSpringel fc Farraii 120071: 
Merten et al.l 120 lit) we calculate their times-since- 
coUision using our analytic model and find that DLSCL 
J0916.2-f 2951 is probing a merger stage > 2 — 5 times 
further progressed than any other known dissociative 
merger (Figure E]). 

6. DARK MATTER IMPLICATIONS 

Given the evident merg er scenario we ar e able to use 
the first method of Marke vitch et al.l (12004) and place 
a rough limit on the DM self-interaction cross-section, 
cdm- This method compares the scattering depth of the 
dark matter, tdm = o'DM'TiQj^SDM, with that of the ICM 
gas, TicM ~ 1, where ttidm is the DM particle mass and 
Edm is the surface mass density of the DM particles. 
Edm is approximately the WL measured surface mass 
densi ty, E, since ^ 80% of a typical cluster's mass is 
DM (iDiaferio et al.ll2008D. For ease o f co mparison with 
the results of IMarkevitch et al.l (|2004f ) and lMerten et alJ 
(|20lll) we examine the surface density averaged over 
the face of the subcluster within r=125kpc, which is 
E fa 0.15 gcm~^; thus we find ctdm'tiqJ^ < 7cm^g~^. 
Note that we cannot apply the velocity-dependent anM 
constraint methods outlined bv IMarkevitch et all ()2004l ) 
since our analytic model assumes ctdm = 0. 

7. DISCUSSION 

While we use DLSCL J0916.2-f 2951 to provide fur- 
ther evidence for the canonical DM model and indepen- 
dently constrain ctdm^tt-dm' believe that its greatest 
value is as a probe for a new and special phase of clus- 
ter formation. It provides a greatly improved temporal 
lever-arm with which to guide numerical simulations that 
explore the major merger phase. This is potentially im- 
portant given that much of our knowledge of the cluster 
merger proce ss comes from num erical hydrodynamic sim- 
ulations (e.g. lPoole et al.ll2006l ). which are used to place 
the tightest c onstraints on ctdm^Idm 0.7cm^g~^; 
iRandall et al.l [20081 and bring observed merger veloci- 
ties (inferred from the observed shoc k velocity) more in 
line with the expectati ons of ACDM ([Springel fc Farrail 
I2007t iLee fc Komatsul [2010). Secondly, the large pro- 
jected separation relative to the virial radii of the sub- 
clusters (-R200 ~ IMpc) enables the deconvolution of the 
subclusters from the Central region and direct compari- 
son of the physical properties of each. This will provide 
new insight into the behavior of the cluster constituents 
(gas, galaxies, & DM) during a major merger. For ex- 
ample, it is well established that galaxy clusters play an 
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Fig. 4.— The phase space of DLSCL J0916.2-I-2951 calculated 
from a Monte Carlo analytical treatment of the system assum- 
ing the clusters to be uniform-density spheres with -R200 radii. 
The measured velocities of the other mergers (with the excep- 
tion of "Bullet SF07"; Bradac ct al. 2008; Mark evitch et al.|[2005l : 
ISpringel & Farrar 2007; Morten ct al. 2011) were input into the 
model to calculate their time-since-collision. It is clear that DLSCL 
J0916. 2-1-2951 probes a significantly different epoch from the ex- 
isting dissociative mergers. The close agreement between our an- 
alytical estimate, "Bullet Analytic" , and the N-body simulation 
ifSpringel fc Farrar.. 2007. 1. "Bullet SF07", lends confidence to the 
results of our analytical treatment. 

important role in the evolution of their member galax- 
ies, but it is still unclear whether cluster merger s trig- 
ger star formation (e.g. [M iller & Owen 2003; Owen et alJ 
12005'; Fer rari et all [20051: IHwang fc Led 12003), quench 
it (P oggianti et al. I2004D . or have no immediate effect 
(Chu ng et al.l 12010 ) . This cluster, with its estimated 
time-since-collision approximately equal to the lifetime of 
A-stars, holds promise for discriminating between slow- 
working processes (e.g. galaxy harassment or strangu- 
lation) and fast-acting process (e.g. ram pressure strip- 
ping). 
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